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Abstract 



Oh; 

An analytical model that describes the interactive buckling of a thin- walled I-section 
strut under pure compression based on variational principles is presented. A formulation 
combining the Rayleigh-Ritz method and continuous displacement functions is used to 
derive a system of differential and integral equilibrium equations for the structural com- 
ponent. Numerical continuation reveals progressive cellular buckling (or snaking) arising 
from the nonlinear interaction between the weakly stable global buckling mode and the 
strongly stable local buckling mode. The resulting behaviour is highly unstable and when 
the model is extended to include geometric imperfections it compares excellently with some 
recently published experiments. 
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The buckling of struts and columns represents the most common type of structural 

instability problem [l|. However, when the compression member is made from slender 
metallic plate elements they are well known to suffer from a variety of different elastic 
instability phenomena. In the current work, the classic problem of a strut under axial 
compression made from a linear elastic material with an open and doubly-symmetric cross- 
section an "I-section" |2|, |3| - is studied in detail using an analytical approach. Under 
this type of loading, long members are primarily susceptible to a global (or overall) mode 
of instability namely Euler buckling, where flexure about the weak axis occurs once the 
theoretical Euler buckling load is reached. However, when the individual plate elements 
of the strut cross-section, namely the flanges and the web, are relatively thin or slender, 
elastic local buckling of these may also occur; if this happens in combination with the 
global instability, then the resulting behaviour is usually far more unstable than when the 
modes are triggered individually. Recent work on the interactive buckling of struts include 
experimental and finite element studies [4J, |5|], where the focus was on the behaviour of 
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struts made from stainless steel. However, the more generic finding that the members had 
an increased sensitivity to imperfections was highlighted. Other structural components 
that are known to suffer from the interaction of local and global instability modes are 
thin- walled beams under uniform bending |6J , sandwich struts [7|, |8J , stringer-stiffened and 
corrugated plates P, LL0| and built-up, compound or reticulated columns (ill ]. 

Apart from the aforementioned work where some numerical modelling was presented |5|] , 
the formulation of a mathematical model accounting for the nonlinear interactive buckling 
behaviour has not been forthcoming. The current work presents the development of a 
variational model that accounts for the mode interaction between global Euler buckling and 
local buckling of a flange such that the perfect and imperfect elastic post-buckling response 
of the strut can be evaluated. A system of nonlinear ordinary differential equations subject 
to integral constraints is derived and is solved using the numerical continuation package 
Auto [12|] . It is indeed found that the system is highly unstable when interactive buckling is 
triggered; snap-backs in the response, showing sequential destabilization and restabilization 
and a progressive spreading of the initial localized buckling mode, are also revealed. This 



or 



latter type of response has become known in the literature as cellular buckling |13f 
snaking |14[ and it is shown to appear naturally in the numerical results of the current 
model. As far as the authors are aware, this is the first time this phenomenon has been 
captured analytically in struts undergoing Euler and local buckling simultaneously. Similar 
behaviour has been discovered in various other mechanical systems such as in the post- 
buckling of cylindrical shells [15J , the sequential folding of geological layers |16J and most 
recently in the lateral buckling of thin- walled beams under pure bending [6J . 

Experimental results from the literature |4j, LL7| are used for validation purposes. The 
mechanical destabilization and the nature of the post-buckling deformation compare ex- 
cellently with the current model. This demonstrates that the fundamental physics of this 
system is captured by the analytical approach both qualitatively and quantitatively. A brief 
discussion is presented on how the current model could be enhanced and then conclusions 
are drawn. 

2. Analytical Model 

Consider a thin- walled I-section strut of length L made from a linear elastic, homoge- 
neous and isotropic material with Young's modulus E and Poisson's ratio v. It is loaded 
by an axial force P that is applied at the centroid of the cross-section, as shown in Figure 
[T](a) and (b) respectively, with rigid end plates that transfer the force uniformly to the 
entire cross-section. The web is assumed to provide a simple support to both flanges and 
not to buckle locally under the axial compression, an assumption that is justified later. 
In the current study, the total cross-section depth is h with each flange having width b 
and thickness t. It is assumed currently that the I-section is effectively made up from two 
channel members connected back-to-back; hence, the assumption is that the web thickness 
t w = 2t, a type of arrangement that has been used in recent experimental studies |4|, [y, \\J_ . 



The strut length L is varied such that in one case, which is presented later, Euler buckling 
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Figure 1: (a) Elevation of an I-section strut of length L that is compressed axially by a force P. The 
lateral and longitudinal coordinates are y and z respectively, (b) Cross-section of strut; the transverse 
coordinate is x. (c) Sway and tilt components of the minor axis global buckling mode, (d) Local buckling 
mode: out-of-plane flange displacement functions Wi(x,z); note the linear distribution in the x direction. 



about the weaker y-axis occurs before any flange buckles locally and in the other case the 
reverse is true - flange local buckling is critical. 

The formulation begins with the definitions for both the global and the local modal 
displacements. Timoshenko beam theory is assumed, meaning that the effect of shear is 
not neglected as in standard Euler-Bernoulli beam theory. Although it turns out that the 
effect of shear is only minor, it is necessary to account for it since it provides the key terms 
within the total potential energy that allow buckling mode interaction to be modelled |6|,lZ|. 
To account for shear, two generalized coordinates q* and q t) defined as the amplitudes of 
the degrees of freedom known as "sway" and "tilt" [7] are introduced to model the global 
mode, as shown in Figure [U(c), where the lateral displacement W and the rotation 6 are 
given by the following expressions: 



W(z) = q s L sin — , 6{z) = q t 7r cos — . 



(i) 



For the present case, the shear strain in the xz plane, j xz , is included and is given by the 

following expression: 

dW „ , 7TZ 



Ixz = — — - = (q s - q t ) 7i cos — . 



(2) 



Of course, Euler-Bernoulli beam theory would imply that since ^y xz = 0, then q s = q t . 

The local mode is modelled with appropriate boundary conditions. Moreover, the 
possibility of a distinct local buckling mode occurring before global buckling implies that 



the entire flange may deflect. However, if the interaction between local and global modes 
occurs then the symmetry of the local buckling mode would be broken and the flanges 
would not buckle with the same displacement. Hence, two separate lateral displacement 
functions w\ and w 2 need to be defined, as shown in Figure [H(d) , to allow for the break in 
symmetry. Since the outstands of the flanges have free edges, whereas the web is assumed 
to provide no more than a simple support to the flanges, a linear distribution is assumed 
in the x direction; Bulson [18| showed this distribution is correct for the local buckling 
eigenmode for that type of rectangular plate. For the local mode in-plane displacements 
u i) the distributions are also assumed to be linear in x, as shown in Figure El This is in 
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Figure 2: Displacement functions of local buckling mode in flanges. Longitudinal and lateral flange 
displacements Ui{x,z) and Wi(x,z) respectively. Note the linear distributions in x direction and the 
average end-displacement, as opposed to the modelled flange end-displacement, which is used to calculate 
the local contribution to the work done. 



fact another consequence of the Timoshenko beam theory assumption where plane sections 
are assumed to remain plane. These assumptions lead to the following expressions for the 
local out-of-plane displacements Wi with the in-plane displacements u{. 
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Wi(z), Ui(x,z) = (-If 
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(3) 



where i = {1,2} throughout the current article. The transverse in-plane displacement 
v(z, x) is assumed to be small and is hence neglected for the current case; this reflects the 
findings from Koiter and Pignataro [9] for rectangular plates with three pinned edges and 
one free edge. 

Since, in practice, perfect geometries do not exist, an initial out-of-straightness in the 
x-direction, Wo, is introduced as a global imperfection to the web and flanges in the current 
model. An initial rotation of the plane section 9 is also introduced to simulate the out-of- 
straightness in the flanges. The expressions for Wq and 9 are given by: 
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(4) 



and are analogous to Equation (pp) . Note that the assumption of Timoshenko beam theory 
implies that shear strains in the xz plane due to the initial imperfection are also introduced. 

2.1. Total potential energy 

The total potential energy, V , was determined with the main contributions being the 
global and local bending energy Ufa and Uu respectively, the membrane energy U m , and 
the work done PS. Note that the global bending energy U\> only comprises the bending 
energy stored in the web, since the membrane energy stored in the flanges accounts for 
the effect of bending in the flanges through the tilt mode. The initial out-of-straightness 
Wq(z) is stress-relieved [113, 120J, implying that the elemental moment M drops to zero as 
illustrated in Figure [3](a). The global bending energy involves the second derivative of W 
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Figure 3: Introduction of geometric imperfections Wo and 9$ in (a) the web and (b) the flanges. 
and Wq and is hence given by: 
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where dots represent differentiation with respect to z and I w — t^(h — 2t)/12 is the second 
moment of area of the web about the global weak axis. Obviously, for the case where 
t w = 2t, the expression becomes I w = 2t 3 (h — 2t)/3. The local bending energy, accounting 



for both flanges, is determined as: 
U bl 
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dz, 



where £>^, the contribution from wi to the standard expression for the incremental strain 
energy from bending a plate [1], is given by: 
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with D = Et 3 / [12(1 — v 2 )} being the plate flexural rigidity. The buckled configuration of 
the flange plate involves double curvature in the x and z directions, indicating the non- 
developable nature of plate deformation. The so-called membrane strain energy (U m ) is 
derived from considering the direct strains (a) and the shear strains (7) in the flanges thus: 
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where: 
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(9) 



The transverse component of strain e xi is neglected since it has been shown that it has no 

effect on the post-buckling behaviour of a long plate with three simply-supported edges 

and one free edge [9] . The longitudinal strain e z has to be modelled separately for different 

outstand flanges. Recall that the tilt component of the in-plane displacement from the 

global mode, including the initial imperfection, is given by u t = —(9 — 9 )x as shown in 

Figure [3](b) ; hence: 

du t 7T 2 . 7TZ 

^global = ~^~ = x (q t - q t0 ) — sin — . (10) 



dz 



The local mode contribution is based on von Karman plate theory. A pure in-plane com- 
pressive strain A is also included. The direct strains in the compression and tension side 
of the flanges, denoted as e z \ and e Z 2 respectively, are given by the general expression: 
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The strain energy from direct strains (Ud) is thus, assuming that /i> t: 
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where, apart from the final term which represents the energy stored in the web, the contri- 
butions are from the direct strains in both flanges. The shear strain energy U s contains the 
shear modulus G ) which is given by E/[2(l + v)] for a homogeneous and isotropic material. 
The shear strain ^y xz contributions are also modelled separately for the compression and 
the tension side of the flanges. The expression for each outstand is given by the general 
expression: 



oz ox oz ox 

TTZ /2\ Ax 

= (Qs ~Qt- QsO + fto) 7TCOS — + (-1)' f - J Ui + -pWiWi. 

The expression for the strain energy from shear is thus: 
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Finally, the work done by the axial load P is given by: 
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(15) 



where £ comprises the longitudinal displacement due to global buckling, the in-plane dis- 
placement due to local buckling and the initial end shortening. Note that the displacement 
due to local buckling is taken as the average value between the maximum in-plane displace- 
ment in the more compressed outstand u\ and the maximum in-plane displacement in the 
less compressed outstand u 2) which is illustrated in Figure [2 Moreover, the possible term 
in q sQ has been neglected since it would vanish on differentiation for equilibrium anyway. 
The total potential energy V is therefore assembled thus: 



V = U bo + U bl + U m - PS. 



(16) 



2.2. Variational Formulation 

The governing differential equations are obtained by performing the calculus of varia- 
tions on the total potential energy V following a well established procedure that has been 
detailed in [7]. The integrand of the total potential energy V can be expressed as the 
Lagrangian (C) of the form: 
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The first variation of V , which is denoted as SV ', is given by: 
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To find the equilibrium states, V must be stationary, which requires SV to vanish for any 
small change in W{ and Ui. By assuming that Swi = d(5wi)/ dz : Swi = d(Swi)/dz and sim- 
ilarly 5iti = d(5ui)/dz, integration by parts allows the development of the Euler-Lagrange 
equations for wi and Ui] these comprise fourth order ordinary differential equations (ODEs) 
in terms of W{ and second order ODEs in terms for U{. For the equations to be solved by the 
continuation package Auto, the system variables need to be rescaled with respect to the 
non-dimensional spatial coordinate z = 2z/L. Non-dimensional out-of-plane displacements 
Wi and in-plane displacements Ui are also introduced as 2wi/L and 2ui/L respectively. Note 
that these scalings assume symmetry about the midspan and the differential equations are 
solved for half the length of the strut; this assumption has been shown to be perfectly 
acceptable for cases where the global buckling is critical [20]. For cases where local buck- 
ling is critical, this condition is also acceptable so long as the length of the strut L is 
much larger than the flange outstand width 6/2; hence the critical loads for symmetric and 
antisymmetric modes are sufficiently close for the buckling plate. The non-dimensional 
differential equations for Wi and Ui are thus: 
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where i = {1,2} again along with D = EtL 2 /D, G = GtL 2 / D and <j) = L/b. Equi- 
librium also requires the minimization of the total potential energy with respect to the 
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generalized coordinates q s , q t and A. This essentially provides three integral conditions, in 
non-dimensional form: 
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where s = 2GtbL 2 /(EI W ) and t = 12G(j) 2 / E. Since the strut is an integral member, the 
expressions in Equation (ED) provide a relationship linking q s and q t before any interactive 
buckling occurs, i.e. when wi = U{ = 0. This relationship is assumed to hold also between 
q s0 and q t0 , which has the beneficial effect of reducing the number of imperfection amplitude 
parameters to one. The relationship between q s0 and q t0 is given by: 

q s o= (j + ljqto- (22) 

The boundary conditions for vbi and Ui and their derivatives are for pinned end condi- 
tions for x = and for symmetry at x = 1: 



Wi(0) = Wi(0) = Wi(l) = Wi(l) = Ui(l) 
with further conditions from matching the in-plane strain: 
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(23) 



(24) 



Linear eigenvalue analysis for the perfect strut (q s0 = q t0 = 0) is conducted to determine 
the critical load for global buckling P C . This is achieved by considering that the Hessian 
matrix V st at the critical load is singular. Hence: 
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Recalling of course that in fundamental equilibrium for this case, q s = q t 
Hence, the critical load for global buckling is: 



(25) 
Wi = Ui = 0. 

- v + TTW (26) 

If the limit G — > oo is taken, which represents a principal assumption in Euler-Bernoulli 
bending theory, the critical load expression converges to the Euler buckling load for an 
I-section strut buckling about the weak axis. 
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3. Numerical examples of perfect behaviour 

The full nonlinear differential equations are obviously complicated to be solved analyt- 
ically. The continuation and bifurcation software Auto-07p J12J has been shown in the 
literature p, [3| to be an ideal tool to solve the equations numerically. For this type of 
mechanical problem, one of its major attributes is that it has the capability to show the 
evolution of the solutions to the equations with parametric changes. The solver is very 
powerful in locating bifurcation points and tracing branching paths as model parameters 
are varied. To demonstrate this, an example set of cross-section and material properties 
are chosen which are shown in Table [TJ In this example, perfect behaviour is assumed and 



Flange width b 96 mm 

Flange thickness t 1.2 mm 

Cross-section depth h 120 mm 

Cross-section area A 513 mm 2 

Young's modulus E 210 kN/mm 2 

Poisson's ratio v 0.3 



Table 1: Cross-section and material properties of an example strut. Recall that the thickness of the web 
t w = 2t. The geometric properties are similar to those tested in |4|. The length L is varied such that the 
cases where global buckling or local buckling are critical can be presented. 

hence W — #o = 0. The global critical load P C can be calculated using Equation f |26l) . 
whereas an estimate for the local buckling critical stress af can be evaluated using the 
well-known plate buckling formula erf = kDn 2 /{b 2 t) ) where the coefficient k depends on 
the boundary conditions; approximate values of k — 0.426 and k — 4 are chosen for the 
rectangular plates representing the flange outstands (three edges pinned and one edge free) 
and the web (all four edges pinned) respectively, assuming that the plates are relatively 
long [18j . Table [2] summarizes the critical stresses and shows that the assigned cross-section 

L (m) a% (N/mm 2 ) ^ c flanse (N/mm 2 ) <^ c web (N/mm 2 ) Critical mode 
3.5 58.3 51.1 2731 Local (flange) 
4.0 447 5L1 2731 Global 

Table 2: Theoretical values of the global and local critical buckling stresses for two separate lengths. The 
expression for a^ = P^/A and the web is obviously not vulnerable to local buckling. 

dimensions satisfy the assumptions that the local mode is critical for one of the lengths and 
the global mode is critical for the other. Moreover, the critical stress of the web is orders 
of magnitude higher than that of the flange, which justifies the assumption stated earlier. 
It should be emphasized that the local buckling critical stress is calculated numerically in 
Auto and is usually marginally higher than the value given in Table El with an error not 
exceeding 5% with the theoretical expressions given above for the long plates. 

Numerical continuation was performed in Auto for the cases where local buckling and 
global buckling were critical in turn. The principal parameters used in the continuation 
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process were interchangeable, but generally q s was varied for computing the equilibrium 
paths for the distinct buckling modes and P was varied for evaluating the interactive buck- 
ling paths. For the case of local buckling being critical, the continuation process initiated 
from zero load with the local buckling critical load Pp being obtained numerically. The 
post-buckling path was then computed by using the branch switching facility within the 
software and the distinct local buckling equilibrium path was computed until a secondary 
bifurcation point S was found. It was from this point that the interactive buckling path 
was found, again through the use of branch switching. For the case where global buckling 
was critical, since the critical load was determined analytically in Equation f |26l ). the initial 
post-buckling path was computed first from Pp and many bifurcation points were detected 
on the weakly stable post-buckling path; the focus being on the one with the lowest value of 
q S) the secondary bifurcation point S. A subsequent run was then necessary starting from 
S using the branch switching function, after which the equilibrium path again exhibits the 
interaction between the global and the local modes. Figure H] shows the procedures for the 
cases diagrammatically with (a-b) concerning the perfect cases discussed above and (c), 
the imperfect case, is considered later in the section on validation. 
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Figure 4: Numerical continuation procedures, (a) Local buckling being critical; (b) Global buckling being 
critical; (c) Imperfect case - example of local buckling being critical shown. The thicker line shows the 
actual solution path in each of the examples shown. Points C and S represent the critical and secondary 
bifurcations respectively, whereas the point So represents the bifurcation leading to interactive buckling in 
the imperfect case with the load Pjj being the ultimate load in the imperfect case. 



3.1. Local buckling critical 

In this section, the strut with properties given in Tabled] with length L being 3.5 m is 
analysed, where the flanges buckle locally first. Figure [5] shows a plot of the normalized 
axial load p = Pj 'Pp versus (a) the global mode and (b) the local mode amplitudes; (c) 
shows the local and global mode relative magnitudes during post-buckling and (d) shows 
that there is a small but importantly, non-zero shear strain during global buckling. The 
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Figure 5: Numerical equilibrium paths for L = 3.5 m where local buckling is critical. Graphs of the normal- 
ized force ratio p versus (a) the generalized coordinate q s and (b) the maximum out-of-plane displacement 
of the buckled flange plate, in non-dimensional form, w max /t are shown, (c) shows w max /t versus q s and 
(d) shows the relationship between the generalized coordinates q s and q t defining the global buckling mode 
during interactive buckling, with the dot-dash line showing the Euler-Bernoulli bending condition q s = q t . 

local critical buckling load is calculated at p — 0.905, whereas according to Table El this 
value should be 0.877, which represents a small error of 3%, particularly since it is well 
known that the theoretical expression for the critical buckling stress for the long plate is 
usually an underestimated value. 

One of the most distinctive features of the equilibrium paths, as shown in Figures E](a)- 
(c), is the sequence of snap-backs that effectively separates the equilibrium path into 10 
individual parts (or cells) in total as shown. The fourth, seventh and the tenth paths are 
labelled as (74, C7 and C10 respectively. Each path or cell corresponds to the formation of 
a new local buckling displacement peak or trough. Figure [6] illustrates the corresponding 
progression of the numerical solutions for the local buckling functions w\ and u\ from cell 
C\ to (77, where C\ represents the initial post-buckling equilibrium path generated from 
C. Once a secondary bifurcation is triggered at S, it is observed that the local buckling 
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Figure 6: Numerical solutions for the local out-of-plane displacement w\ (left) and local in-plane displace- 
ment u\ (right) for the tip (x = —b/2) of the vulnerable flange. Individual solutions on equilibrium paths 
CI to CI are shown in sequence from top to bottom respectively. 

mode is contaminated by the global mode and interactive buckling ensues with the buckling 
deformation spreading towards the supports as new peaks and troughs are formed. Figure 
[3 shows a selection of 3-dimensional representations of the deflected strut that comprise 
the components of global buckling (W and 6) and local buckling (w{ and Uj) at a specific 
state on paths CI, C4, C7 and C10. As the equilibrium path develops to C10, the 
maximum out-of-plane displacement w max approaches a value of 2.5 mm which is roughly 
twice the flange thickness and can be regarded as large in terms of geometric assumptions. 
The interactive buckling pattern becomes effectively periodic on path C10. Any further 
deformation along the equilibrium path would be expected to cause restabilization to the 
system since the boundaries would begin to confine the spread of the buckling deformation. 
It should be stressed of course that any plastic deformation during the loading stage would 
destabilize the system significantly. Figure [8] shows the comparison between the lateral 
displacement of the two flange outstands. The local buckling displacement in the non- 
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Figure 7: Numerical solutions of the system of equilibrium equations visualized on 3-dimensional repre- 
sentations of the strut. The results are shown for individual points on paths (a) CI (p = 0.9039), (b) 
C4 (p = 0.9081), (c) CI (p = 0.8111) and (d) C10 (p = 0.7177). All dimensions are in millimetres. 
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Figure 8: Numerical solutions for the local out-of-plane displacement w 2 (left) for the tip of the non- 
vulnerable flange (x = b/2) and w\ (right) for the tip of the vulnerable flange (x = —b/2) for cells 1-3. 
Note the rapid decay of w 2 reflecting the reducing compression in that outstand once global buckling is 
triggered. 

vulnerable outstand w 2 decays to zero rapidly as the global mode amplitude increases 
during interactive buckling; by the third cell, w 2 has vanished implying that if global 
buckling occurs first, both w 2 and u 2 would be negligible. 

The magnitude of direct and shear strains may be calculated once the governing dif- 
ferential equations are solved. The direct strain in the non- vulnerable part of the flange 
becomes tensile at (710 due to bending, whereas the maximum direct strain in the vulnera- 
ble part of the flange is approximately 1.3 x 10 _3 (= 0.13%). This level of strain is confined 
to the ends of the strut and is also well below the yield strain of most structural steels; 
moreover for the stainless steels given in Becque and Rasmussen |4j, LL7J , significant strain 
softening only begins from approximately 0.15% strain and so quantitative comparisons 
can be made for the post-buckling response for the majority of the cells. 

Systems that exhibit the phenomenon described above are termed in the literature to 
show "cellular buckling" |13|] or "snaking" [14]. In such systems, progressive destabiliza- 
tion and restabilization is exhibited; currently, the destabilization is caused primarily by 
the interaction of the global and local instabilities, whereas the restabilization is caused 
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by the stretching of the buckled plates when they bend into double curvature. As the 
amplitude of the global buckling mode q s increases, the compressive bending stress in the 
flange outstands increase also, which imply that progressively longer parts of the flange are 
susceptible to local buckling. Since local buckling is inherently stable, the drop in the load 
from the unstable mode interaction is limited due to the stretching of the plate when it 
buckles into progressively smaller wavelengths. Therefore, the cellular buckling occurs due 
to the complementary effects of the unstable mode interaction and stable local buckling. 

3.2. Global buckling critical 

The strut with properties given in Table [U with length L being 4.0 m is now analysed; 
in this case, global Euler buckling occurs first. Figure [9] shows plots of the equilibrium 
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Figure 9: Numerical equilibrium paths for L = 4.0 m where global buckling is critical. Graphs of the 
normalized force ratio p versus (a) the generalized coordinate q s and (b) the maximum out-of-plane dis- 
placement of the buckled flange plate, in non-dimensional form w; max /t, are shown, (c) shows w; max /t 
versus q s and (d) shows the relationship between the generalized coordinates q s and q t defining the global 
buckling mode during interactive buckling, with the dot-dash line showing the Euler-Bernoulli bending 
condition q s = q t . 

diagrams that correspond directly to Figure El This time, cellular buckling is triggered 
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when the pure global mode is contaminated by the local mode. Since the global mode is 
only weakly stable, no significant post-buckling stiffness is exhibited initially. Moreover, 
since the global mode places the non- vulnerable flange outstand into less compression before 
any local buckling occurs, the functions w 2 and u 2 can be neglected as a consequence of the 
observations made in connection with Figure [8j this simplifies the formulation considerably. 
The emergence of the buckling cells in sequence is very similar to that shown for the 
case where local buckling is critical and so it is not presented in detail for brevity. Nev- 
ertheless, with the model in place, quantitative comparisons can be made against existing 
experiments. 

4. Validation and discussion 

4-1- Comparison with experiments of Becque and Rasmussen 

A recent experimental study of thin- walled I-section struts by Becque and Rasmussen 
[4J, lilt] focused on the case where local buckling is critical. Although the struts were made 
from a stainless steel alloy (ferritic AISI404), the compressive stress-strain curve showed 
that the material remained linearly elastic when the strain was below approximately 0.15%. 
Two specific tests were conducted on struts with material and geometric properties as given 
in Table [31 The initial out-of-straightness mid-length lateral deflections of the specimens of 



Strut length L 


3.0 m 


2.5 m 


Flange width b 


96.64 mm 


96.80 mm 


Corner Radius r 


3.06 mm 


3.02 mm 


Flange thickness t 


1.21 mm 


1.21 mm 


Section depth h 


125.12 mm 


125.24 mm 



Table 3: Geometric properties for the strut tests taken directly from [4J, LL7j. Recall that the thickness 
of the web t w = 2t. For both struts the initial Young's modulus E = 195 kN/mm 2 and Poisson's ratio 
v = 0.3. 

length L being 3 m and 2.5 m were measured to be L/3352 and L/16234 respectively [171 ]. 
In order to make direct comparisons, numerical runs were conducted in Auto with the 
initial global buckling mode imperfection amplitude ratio q s0 being equal to 3 x 10 -4 and 
6 x 10 -5 respectively. The cross-section properties given in Table [3] were adapted slightly 
to consider the effective width of the flange, b e = b — 2r; the effective width was used in 
the numerics for the analytical model, the results of which follow. 

The numerical continuation process was initiated from zero load with the process being 
illustrated in Figure H](c). The value of q s was increased up to a bifurcation point, shown 
as So in Figure [4](c), after which interactive buckling was introduced. The equilibrium 
path then progressed to a limit point at which P can be defined as the ultimate load 
Pjj. Then destabilization and the cellular buckling behaviour was observed as described in 
the previous section. Figures [10] and | 
model, the experimental results from 
comparisons show strong agreement between the analytical model and the results from 

17 



[TT show comparisons between the current analytical 
l [4J, llTI] and the numerical models from [5|, LL7|] . The 
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Figure 10: Numerical equilibrium paths comparing with Becque's experiment with a strut of length 3 m. 
Graphs of the applied axial load P versus (a) the total end shortening, (b) the generalized coordinate 
q s , (c) the generalized coordinate q t , (d) the out-of-plane displacement of the buckled flange plate w\ 
measured at z — 400 mm are shown. Solid lines show the current analytical model, whereas the dashed 
and dot-dashed lines respectively show the experimental and finite element results from [17j. 

the physical experiments, the correlation being clearly superior to the previous numerical 
results presented in p, [L7J . 

For the 3 m length strut, the ultimate load was found to be 25.2 kN from the experiment, 
which is approximately 3% higher than the numerical value from the analytical model where 
Py = 24.4 kN. It is also observed in Figure [TUlfd) that the theoretical local out-of-plane 
displacement w\, at a location that was remote from the strut midspan (z = 400 mm), 
changes from positive to negative and vice versa several times. This is a signature of the 
cellular behaviour, indicating the progressive change in wavelength of the local buckling 
mode pattern. The actual experimental response, on the other hand, is always likely to 
jump to the final cell relatively rapidly once the initial instability is triggered. This was 
shown in the experiments conducted during work on the interactive buckling of beams 
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Figure 11: Numerical equilibrium paths comparing with Becque's experiment with a strut of length 2.5 m. 
Graphs of the applied axial load P versus (a) the total end shortening, (b) the generalized coordinate 
g s , (c) the generalized coordinate q t , (d) the out-of-plane displacement of the buckled flange plate w\ 
measured at z — 400 mm are shown. Solid lines show the current analytical model, whereas the dashed 
and dot-dashed lines respectively show the experimental and finite element results from [17j. 



, particularly in the cases where global and local buckling were triggered at similar load 
levels. The reason for this is that in a physical experiment, even with displacement control, 
the mechanical response in the region with snap-backs exhibits dynamic, rather than static, 
behaviour. Although in the current case the experiment did not pick up the full cellular 
response, it did show the change from positive to negative for w\, which is a clear indication 
of the changing wavelength in the local buckling mode pattern. The interactive buckling 
wavelength A can also be compared, which is defined in Figure [121 The local buckling 
mode had a plate buckling wavelength that was measured to be 275 mm with a modulated 
amplitude for this specific test [4] . The numerical results in the current work show that the 
value of A is 280 mm for the interactive buckling wavelength at the end of the equilibrium 
paths from the analytical model shown in Figure [101 The close comparison (less than 2% 
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Figure 12: Definition of local buckling wavelength A from results for w(= w\) from the variational model. 

difference) offers further grounds for encouragement for future developments of the current 
model. 

For the 2.5 m length strut, the features are similar; the ultimate load Pjj is 16% higher 
than the maximum load shown in the experiment. However, this is only a very small 
part of the global picture. The graphs in Figure ITTlf a-c) show sequential snap-backs in 
the theoretical response almost immediately in the post-buckling range that reduce the 
true load-carrying capacity to levels which practically coincide with the experimental re- 
sult, demonstrating an excellent overall comparison. In Figure [TT](d), a similar response is 
observed to Figure [TUlfd). The change in sign of the out-of-plane displacement from the 
test clearly demonstrates the changing buckling wavelength again. Unfortunately, a nu- 
merical measurement of the plate buckling wavelength was not reported for this particular 
experiment. 

Finally, it can be seen that the results from the analytical model and the experiments 
begin to diverge after a certain level of displacement. This is postulated to be as a result 
of material softening due to the use of stainless steel in the experiments, whereas linear 
elasticity is assumed throughout the analytical model. However, for the most part the close 
comparisons between the analytical model and the experimental results, in the authors' 
opinion, validate the current modelling approach both qualitatively and quantitatively. 

4-2. Future model enhancements 

The success in capturing the interactive buckling behaviour allows for some speculation 
of how the current work may be extended. The technical difficulty of capturing sharp 
successive snap-back instabilities numerically most probably explains why the previous 
finite element models |5|, although giving safe predictions for the global strength of the 
tested struts in jj, ll7|, showed a relatively indifferent comparison with the experiments. A 
possible way around this problem in future finite element modelling of such struts could be 
to introduce in turn a sequence of initial imperfections with different shapes that resemble 
the modes from each cell. This would allow an envelope of the nonlinear equilibrium 
solutions to be computed that resemble the actual post-buckling response. 

Other issues that can be investigated are those regarding the assumption of the fixity 
between the elements of the cross-section, namely between the web and the flanges, and the 
introduction of lipped ends to the flanges. In terms of joint fixity, the flange-web junctions 
are modelled as pinned and hence are free to rotate. By modelling them as partially to 
fully rigid, a more extensive range of responses could be captured. With the resultant 
increase in structural stiffness in the cross-section this would introduce to the system, the 
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local buckling load would definitely increase. However, the early evidence from a pilot 
study is that the post-buckling becomes less cellular as a result [21]. A similar effect may 
be obtained by attaching or designing flanges with lips to reduce the vulnerability to local 
buckling. However, lips introduce the possibility of distortional buckling [22] which, in this 



context, is known to resemble localized buckling [23l. 1241] rather than the cellular buckling 
found presently. Localization, in this context, can be more severely destabilizing than the 
cellular buckling. Work on this latter enhancement is currently in the early stages and 
hence the point regarding the potentially greater severity in the post-buckling instability 
is purely conjecture currently. 

5. Concluding remarks 

A nonlinear analytical model based on variational principles has been presented for 
axially-loaded thin-walled I-section struts buckling about the weak axis of bending. The 
model identifies an important and potentially dangerous interaction between global and 
local modes of instability, which leads to highly unstable cellular buckling through a series 
of snap-back instabilities. These result from the increasing contributions of buckling mode 
amplitudes forcing the flanges in more compression to buckle progressively. This process 
had also been observed in recent experimental work and in other components that suffer 
from a nonlinear interaction between global and local buckling. Comparisons with pub- 
lished experiments are excellent and validate the model. Extending the analytical approach 
would allow further study of the parameters that drive the behaviour and provide a greater 
and more profound understanding of the underlying phenomena. This, in turn, would pro- 
vide designers with the information about the sensitivity of thin-walled components to 
small changes in geometry. 
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